Two-band Effects in the Angular Dependence of Hc2 of MgB2 Single Crystals 



A. Rydh0 U. Welp, A. E. Koshelev, W. K. Kwok, and G. W. Crabtree 
Materials Science Diviston, Argonne National Laboratory, 9700 S. Cass Ave., Argonne, IL 60439, USA 

R. Brusctti, L. Lyard, and T. Klein 
Lahoratoire d'Etudes des Proprietes Electroniques des Sohdes, CNRS, BP 166, 38042 Grenoble, France 

C. Marcenat 

Departement de Recherche Fondamentale sur la Matiere Condensee, SPSMS, CEA-Grenoble, 38054 Grenoble, France 

B. Kang, K. H. Kim, K. H. P. Kim, H.-S. Lee, and S.-I. Lee 
NCRICS and Dept. of Physics, Pohang University of Science and Technology, Pohang 790-784, Republic of Korea 

(Dated: February 2, 2008) 

The angular dependence of the upper critical field Hc2 of MgB2 single crystals is studied at var- 
ious temperatures by means of specific heat and transport measurements in magnetic fields up to 
17 T. Clear deviations from Ginzburg-Landau behavior are observed at all temperatures and are 
explained by two-band effects. The angular- and temperature dependence of the deviations are in 
qualitative agreement with theoretical predictions based on band-structure calculations. Quantita- 
tive agreement is obtained with an interband coupling slightly stronger than the calculated, enabling 
band-structure anisotropies and interband coupling strength to be experimentally estimated. This 
provides a new pathway to the study of disorder and doping effects in MgB2 . 

PACS numbers: 74.25.Bt, 74.25.Dw, 74.25.Fy, 74.25.0p 



The emergence of new theoretical works with close 
experimental connections has significantly deepened the 
understanding of the properties of magnesium diboride 
(MgB2). Despite the fact that the superconducting 
properties of MgB2 with its fairly simple atomic struc- 
ture were just recently discovered jll, this phonon me- 
diated s-wave superconductor has already been the sub- 
ject of intense and numerous studies due to its exotic 
properties arising from a complex, disconnected, multi- 
band Fermi surface. Band-structure calculations have 
demonstrated that the Fermi surface is composed of pairs 
of three-dimensional 7r-bands and quasi-2D cr-bands 
This effective two-band structure has been confirmed by 
de Haas-van Alphen measurements Q and angle-resolved 
photoemission spectroscopy [^. 

The superconducting properties of the two sets of 
bands are quite different, due to the low overlap of the 
orthogonal cr- and 7r-band wave functions. The supercon- 
ducting gap ranges from 1.5 to 3.5 meV on the 7r-bands 
and from 5.5 to 8 meV on the strongly superconducting 
cr-bands ^ . This double-gap nature has been verified by 
tunneling experiments 7, £.1 , heat capacity measurements 
, and spectroscopy O Q • 

Theoretically, two-band superconductivity has a his- 
tory starting well before MgBz HI Q Through 
theoretical advances, a fairly unified picture has emerged 
with predictions that can be experimentally substanti- 
ated [illllllllllllliililillillillil. One of 
the salient predictions associated with a pronounced two- 
band effect is a difference between the coherence length 
anisotropy 75 = £,ab/£,c 0, 0, 1^ and the pene- 



tration depth anisotropy jx — Xc/Xab HMIlll, both of 
which become temperature dependent with opposite ten- 
dencies. For MgB2, a strong decrease of 75 = H'^2lHc2 
from T^fO) ~ 5 to ifiTr) w 2 is found experimen- 
tally IM IH El S IHO, EH, while controversy remains 
about the experimental temperature dependence of 7a 

ESElEi 

In this Letter we present evidence of clear deviations 
of the angular dependence of Hc2 from the anisotropic 
Ginzburg-Landau (GL) description. The Hc2 {T, 9) tran- 
sition of MgB2 single crystals was determined from re- 
sistivity measurements and specific heat with excellent 
agreement between the two. With a slight adjustment of 
some of the parameters supplied by band-structure cal- 
culations, good quantitative agreement is found between 
theory 21] and experiment, yielding fundamental esti- 
mates of band-structure anisotropies and the interband 
coupling strength. 

Several MgB2 crystals with typical dimensions 50 - 
250 /im were obtained through a high pressure heat treat- 
ment of a mixture of Mg and B in excess Mg as de- 
scribed elsewhere [s^. The crystals had Tc values of 
34 - 36 K and a Hc2(fi) ~ 3.5 T. Transport measure- 
ments were performed using standard AC techniques at 
23 Hz with a resolution better than 0.5 nV. For specific 
heat measurements, the crystals were mounted on top 
of flattened 12.7 fim chromel/constantan thermocouple 
junctions. Small temperature oscillations of the sample 
were induced by either a resistive heater wire (Sample 1, 
Argonne) or by modulating the temperature of the cop- 
per base (Sample 2, Grenoble). The temperature 
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FIG. 1: Top: Resistive Hc2 transition at T = 27.5 K as a 
function of magnetic field H and angle 9 from the basal plane. 
Bottom: Corresponding specific heat signature at T = 25.0 K. 

oscillation was measured through the AC voltage across 
the junction, the absolute base temperature monitored 
by a Cernox thermometer, and the sample temperature 
offset obtained from the DC voltage of the thermocouple. 

Figure Q shows the transitions from resistivity (top) 
and specific heat (bottom) as a function of angle at 
T = 27.5 K and T = 25.0 K, respectively. The resistive 
transitions were measured at a relatively high current 
density to suppress the effects of surface superconduc- 
tivity at the well-shaped crystal surfaces, as discussed 
in [23- The value of Hc2{6) was determined through a 
linear extrapolation of the steep drop to zero resistiv- 
ity as shown by the nearly vertical dashed lines in the 
top panel. The appearance of the peak effect just be- 
low Hc2 for some angles is evident in the figure. The 
thermodynamic signature of Hc2 was defined from the 
midpoint of the specific heat transitions, as illustrated in 
the bottom panel. The choice of definition was checked 
not to be significant. It is interesting to note that the 
specific heat step height is fairly independent of the field 
direction. This is in agreement with GL theory, where 
the step height should scale with T{dHc/dT)'^ , where He 
is the (isotropic) thermodynamic critical field. Possible 



FIG. 2: Angular dependence of the upper critical field at 
12.1 K and 25.0 K. The solid lines correspond to the GL 
theory. Small, but clear and consistent deviations from the 
anisotropic, effective-mass description are seen. 

deviations from a constant step height arising from two- 
band effects are too small to be resolved in the current 
data due to uncertainties in the experimental method. 

Clear deviations from an anisotropic GL description 
are, however, seen in the angular dependence of Hc2- In 
Fig. 13 the Hc2 (0) curves are shown for two selected tem- 
peratures together with corresponding fits to the effective 
mass description H^2^{0) = H^^/^cos^ 9 + -/Isin'^ 9)'^/'^. 
The relative deviations are fairly small at 12.1 K as com- 
pared to 25.0 K. They are nevertheless clearly discernible 
at all temperatures and are reproducible between differ- 
ent measuring setups, samples, and methods. Resistive 
measurements by Eltsev et al. displayed similar devia- 
tions but were not analyzed in detail "36] . Deviations 
were also reported at 33 K using torc^uc measurements 
|37| . The latter, however, suffer from the inability to 
measure Hc2 along the symmetry axes. 

The deviations from GL behavior are illuminated by 
plotting the ratio A = [H c2{9) / Hf^ {9)]^ as a function 
of cos^ 9 as shown in Fig |3| When the field is directed 
along the c axis or within the basal plane there are no de- 
viations, since the experimental ^■nd were used 
as parameters for the GL-fit (.4 = 1) at each tempera- 
ture. The shape of the deviations as a function of angle 
is similar for all temperatures, with a maximum ampli- 
tude at around 9 = 20°-30°, i.e., for cos^ - 0.9. The 
theoretical curves are discussed below. 

To investigate the temperature dependence of the de- 
viations in more detail, the maximum amplitude of A{9) 
is plotted as a function of temperature in Fig. ^ (top). 
Good agreement is found between the transport and spe- 
cific heat data from Argonne (Sample 1, measured on 
the same crystal, Tc = 35.5 K) and specific heat mea- 
surements taken at Grenoble (Sample 2, Tc ~ 34 K), 
illustrating the fundamental and consistent nature of the 
deviations. The amplitude is relatively small at low tem- 
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FIG. 3: Ratio A = [Hc2{e)/ H^2^{e)f as a function of cos^ 61 
for the two temperatures of Fig. |21 The dashed line A — 1 
corresponds to the GL description, assuming an anisotropy 
parameter 7^ = H^2/Hc2 that is allowed to vary as a function 
of temperature. Solid curves are given by two-band theory 
with best-fit parameters as discussed in the text. 



peratures and reaches a maximum slightly below Tc- By 
comparing the temperature dependence of Amax with 
that of 7^ (bottom panel) one can see that the maximum 
of AmaxiT) occurs at intermediate values of 7^. 

The general experimental features of Hc2 {T, 9) can be 
excellently described by the recent theory of the angular 
dependence of dirty two-band superconductors |^] . The 
two-band theory requires as input (i) the matrix of effec- 
tive coupling constants A^^ = X^p — where X^p are 
the electron-phonon coupling constants and /u* ^ are the 
Coulomb pseudopotentials (a and /3 are indices for the 
a- and 7r-bands), (ii) band anisotropics and 77r, and 
(iii) a ratio of the diffusion constants in the two bands, 
e.g., rz — V-^^z/'Da^z- The theoretical dependencies of 
A(T) and 7{(T) obtained by using coupling constants 
and anisotropics supplied by band-structure calculations 
[17I are illustrated by the dashed curves in Fig. 21 
It is clear that, while qualitatively similar, the theoret- 
ical curves are displaced closer to Tc and the predicted 
anisotropy is higher than the experimental one. 

The shapes of the theoretical curves are sensitive 
to mainly two parameters, the interband coupling 
strength, expressed through the reduced parameter 5*12 ~ 
Acr7rA^o-/(Ao-cr — A^tt)^ , and the ratio r^. In partic- 
ular, the overall change of anisotropy can be estimated 
as 7^(rc)/7^(0) ~ \/\/\ + Surz- We found that a quan- 
titative description of the experimental data requires (i) 
a significantly increased relative interband coupling S12, 
(ii) a somewhat decreased anisotropy 7^ of the a-band, 
and (iii) an almost isotropic 7r-band. 6*12 was increased 
by augmenting the two off-diagonal coupling constants 
by a factor 1.9 from the values provided in Ref. [23 |. 
The set of parameters that gives the best description of 
the experimental data is listed in Table together with 
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FIG. 4: Top: Maximum deviations of A{9) from the GL the- 
ory as a function of reduced temperature. The dashed curve 
is taken from |2lll . Bottom: Temperature dependence of the 
experimental anisotropy 7^. Solid circles are taken from from 
[23. Other symbols are as above. The inset shows ^max as a 
function of anisotropy, illustrating maximum deviations at in- 
termediate values of 75 where both bands contribute equally. 
Parameters for the theoretical curves are given in Tabled] 



the original parameters. The resulting fits are shown 
as solid curves in Figs. 13 and ^ We note that the 
fits in Fig. 0] allow for an independent determination 
of V^^^/V^^x = {'^^l^^f Itz ~ 0.23, in good agreeiii«it 
with the observation of an enlarged vortex core |23, l2a| ■ 
A possible source of the discrepancy between the the- 
oretical calculations and the experiments is a theoretical 
overestimation of the off-diagonal Coulomb pseudopoten- 
tials (see discussion in Ref. 39] ) resulting in too low val- 
ues of Acttt and Atto-. Unfortunately, no direct exper- 
imental probe of the off-diagonal coupling constants is 
available at present. On the other hand, the Tc values 
of our single crystals are somewhat lower than for poly- 
crystalline samples. Thus, it is also possible that the 
discrepancy is due to a slightly modified band-structure 
arising from a non-stoichiometric composition of the crys- 
tals. Another origin of Tc suppression in single crystals 
could be interband impurity scattering. A natural ques- 
tion is how interband scattering affects the anisotropic 
properties. By analyzing the theoretical corrections to 
the components of the upper critical fields from the in- 
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TABLE L Parameters used in the theoretical computations. 
Parameter Predicted Params." Best Experimental Fit 
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1.02 




6.3 


5.4 




0.034 
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1^ 0.091 0.278 J 


1^ 0.171 0.278 




300 


120 



"As pre dicte d by band-structure calculations. A -val ues are taken 
from Ref. |24| and 7- values are obtained from Ref. Il7| . The param- 
eter = 300 was not calculated b ut e stimated, and corresponds 

to Va,^/V.^,^ = i^ahn^f/rz = 0.2 



elusion of weak interband seattering, we conclude that it 
is unlikely that this scattering is responsible for discrep- 
ancies between the calculated band-structure and exper- 
iment (e.g., for the lower value of 70-) 

MgB2 single crystals are usually described as fairly 
clean, with the cr-band probably in the clean limit 
0, [33, 113 ■ This is supported by the low, reproducible 
value of 7?^2(0) ~ 3.5 T. The temperature dependence of 
7^ has also been described successfully within the clean- 
limit formalism |2oj . However, this still requires similar 
modifications of coupling and band anisotropics from the 
predicted values. To our knowledge, a clean-limit cal- 
culation of the angular dependence of Hc2 has not yet 
been presented. One could expect deviations from the 
anisotropic GL dependence of Hc2 at low temperatures 
even for a clean cr-band due to Fermi surface effects. On 
the other hand, the theory 21] for Hc2{d) should remain 
valid even for a clean cr-band in the vicinity of , where 
this band is described by GL theory and only the contri- 
bution from the 7r-band requires a microscopic approach. 

In summary, we have studied the angular and tem- 
perature dependence of the upper critical field of MgB2 
single crystals by means of heat capacity (specific heat) 
and transport measurements. Clear two-band effects are 
found in both Hc2{d) and the temperature dependence 
of the upper critical field anisotropy 7f (T) . The exper- 
iments are well explained by the theory 121', providing 
a deep understanding of the microscopic parameters de- 
scribing the system. This work thus points out a new 
pathway to the study of disorder and doping effects in 
MgB2, with great implications for future applications. 
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perature has to be significantly smaller than the relative 



change of Tc- In addition, the correction to the basal- 
plane upper critical field has an extra small factor 1 / y/rZ. 



